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Abstract 

A new Darboux transformation is presented for the Hirota-Satsuma coupled 
KdV system. It is shown that this Darboux transformation can be constructed 
by means of two methods: Painleve analysis and reduction of a binary Darboux 
transformation. By iteration of the Darboux transformation, the Grammian type 
solutions are found for the coupled KdV system. 
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1 Introduction 



Hirota and Satsuma [p] proposed the very first coupled KdV system, which describes 
interactions of two long waves with different dispersion relations. They then constructed 
the three-soliton solutions and five conserved quantities for this system. In a following 
paper 0, these authors shown that the coupled KdV system is the four- reduction of the 
celebrated KP hierarchy and its soliton solutions can be derived from those of the KP 
equation. 

It is known that a solitonic equation normally possesses a Lax pair. Using the 
Wahlquist-Estabrook prolongation procedure, Dodd and Fordy || found a Lax repre- 
sentation for the Hirota-Satsuma coupled KdV (HS-KdV) system. Meanwhile, Wilson |4[] 
observed that the HS-KdV system is just an example of many integrable systems arose 
from the Drinfeld-Sokolov theory. A Backlund transformation is constructed for this sys- 
tem by Levi |J. Very recently, some new solutions for HS-KdV system are presented 

Most integrable equations have Darboux Transformations (DT), which are very effec- 
tive to construct solutions. We notice that Leble and Ustinov J7| found a DT for the 
HS-KdV system. They started with an elementary DT for a more general spectral prob- 
lem, then they found that a proper reduction led to a DT for HS-KdV system. In general, 
an integrable system may possess more than one DT. For example, the KdV equation 
has classical DT and binary DT as well. We are interested in finding a new DT for the 
HS-KdV system. We will show that this new DT appears in two ways, from Painleve 
analysis and from reduction of the binary DT for a more general spectral problem (we 
refer to [§]-[|TI| and the references there for Painleve analysis approach to DT). 



The paper is set out as follows. In section 2 we construct a new DT for the HS-KdV 
system within the framework of Painleve analysis. In section 3, we show that a DT can be 
obtained from the reduction of a general DT, then we prove that this DT is equivalent to 
the one found in section 2. In section 4 we discuss the iteration of the DT and show that 
the solutions for the HS-KdV system can be represented as Grammians. Final section 
presents our conclusion. 



2 DT from Painleve analysis 

The HS-KdV system reads 

u t = -u xxx + 3uu x - 6vv x , v t = -v xxx - 3uv x . (1) 

where the subscripts denote partial derivatives. 

A detailed Painleve analysis of this system was performed by Weiss [ I 1 j [12|. He found 
that this system has the Painleve property. Indeed, there are two branches and one of 
them, principal branch, is 



oc 



U = T 2 WjT 7 , V = T 1 V j T ^ 1 (2) 
j=0 j=0 



and the resonances are 

j = -1,0,1,4,5,6 
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By truncating the expansions (|2|) on the constant level, one has a transformation 



v 



u = u 2 + 2(lnr) xx , v = vi H 



T 



(3) 



where u 2 and v± solve the HS-KdV system ([I]) and r, u 2 , vo, v\ have to satisfy the following 
overdetermined system 

r t + r xxx + 3t x u 2 = 2r x tf, 
v = T X H, 



T x 



H t = - 



i r 3 



H xx + ±H 3 + M + ^{t;x}H 



where {r; x} = (jf 31 ) — § is the Schwarzian derivative. 

To obtain a Lax pair for the HS-KdV system ([[]), Weiss introduced a new variable 



W 



T, 



xx 

r x 



(4) 



and derived the equations for H and W, which are modified Hirota-Satsuma system. The 
Miura map between HS-KdV system and its modification is 



u 2 

Vi 



W x + -W 2 + -H 2 - -v> 



H x + WH+-(\ 2 + $ 2 ) 1 / 2 
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Taking 



W + H = 2 



Vlx 



W-H 



,V2x 
' V2 



■# = Asinho!, a = ln( — ), 

V2 



one arrives at the linear problem 



(5) 
(6) 



Vlxx + («2 + Vt)T]i = -\r}2, 

V2xx + {u 2 - v x )r]2 = Xrji, 



(7) 
(8) 



and 



Vit = ~-( u 2x-2v lx )r] 1 + (u 2 -2v 1 )r] lx -2\r] 2x , 
V2t = ~\{u 2x + 2v lx )r] 2 + (u 2 + 2v 1 )rj 2x + 2\r) lx , 



(9) 
(10) 



where the parameter A is scaled for convenience. This is the linear problem found by 
Dodd and Fordy via the prolongation approach ||. 
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We intend to construct a DT for the HS-KdV systems ([!]). As a first step, we rewrite 
the transformation ([3]) between two solutions of the system (H) in terms of the solutions 
of linear problem (|T|-|T0|). Using (H) and (B), we obtain 



u — U2 = 2(lnr) 



V — V\ 



V2Vlx - V2xVl 



11) 



where (771,772) is a solution of the linear system (|7j-|T0|). The x-derivative of function r is 
defined in terms of rji and rj 2 as follows 



VlV2- 



(12) 



The transformation (|TTj) is the Darboux transformation on the potential level. Next, we 
need to find the transformations for the wave function. Let us suppose that (0i, 02) solves 
the linear systems (|7|-|roD with A = \x. For the new wave function, we make the following 
Ansatz 

01 = 01 + -, 02 = 02 + -, (13) 
T T 

where / = f(x,t) and g = g(x,t) are the functions to be determined. By requiring the 
new wave function (0i,02) solves the linear problem (|7|-|T0|) with replacement u-i — > u, 
V\ — > v and A — > fi, we are led to two equations 



+ («2 + Vi)4>l + /i02 + 
1 



+ - 

T 



fxx + (u 2 + v{)f + 2r^.0i + (f]i x r] 2 - ViV2x)(pi + ng 
1 



+ 



2fxT x + fTxx ~ 2r x 0i + (r] lx 7] 2 - T]lV2x)f 



0. 



(14) 



and 



^2xx + («2 - ^l)02 - yU01 + 
1 



+ - 

T 



9xx + ("2 - ^l)5 + 2r x;c 02 - (771^772 - VlV2x)4>2 - A*/ 



1 r 

+ ^T 



2«/aTa. + ^t-kb - 2r 2 2 - (r] lx r]2 - r)iV2x)g 



0. 



(15) 



Let us consider the equation (|14]) first. The first part of this equation, the coefficient of 
r°, is identically zero. Equating the coefficients of r _1 and r -2 to zero respectively, we 
have 



fxx + (u 2 + vi)f + 2r x:c 0i + (771772 - ViV2x)(f>i + yug = 0, 

-2/^ + /Taa, - 2^01 + (77^.772 - Tjir}2x)f = 0. 



From (|T7|) , we obtain 
Now, using (|IED and ( 



/* - /— + ?7l7720l = 0. 

we have 

Ar? 2 / - Wig = rjmifaVix ~ <PixVi)- 



(16) 
(17) 

(18) 
(19) 
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A similar consideration for the other equation ( |15| ) provides us 

In this way, we obtain a system of linear equations for / and g, namely (II 
Solving this linear system, we get the expressions for / and g as follows 



(20) 
and (M\. 



f 
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A 2 — [i 2 



V((f>lVlx - 4>lxVl) + KV2(p2x ~ V2x4>2) 



(21) 
(22) 



Thus, we obtain the explicit expressions for the transformations for the wave functions, 
that is, the equations (0) with / and g given by (pl] - p2]) . However, this is not the end of 
story: we have to make sure that the transformed wave functions and fields also satisfy 
the temporal part of the linear problem. But first we need to obtain the time evolution 
of r, which is 

r t = 2X(r]f - r]\) - 2r] lx r] 2x - u 2 ri 1 r] 2 . (23) 

Verifying that the new wave functions and new fields obtained under the transformation 
are indeed the solutions for transformed linear systems is too tedious to perform by pen 
and paper, but we do check its validity by means of MAPLE. Thus, we obtain a new DT 
for the HS-KdV system ([!]). We summarize our results as 

Proposition Let (7/1,772) be the solutions of the linear systems (ffll^j and (4>i,<f)2) be 
the solutions of the same linear systems with the spectral parameter /i. Let the function 
t be defined by equations ^U\ ) and i\2^). Then the new field variables and the new wave 
functions defined via (11) and (|73|j with f and g given by (21) and (2/k) solve the linear 
system (0-[T^) with as the spectral parameter. 



3 DT from reduction 

Our DT constructed above is different from the known DT obtained by Leble and Ustinov 
0. Indeed, their DT is of kind classical DT while our one is of kind binary DT: we used 
both the wave functions and the adjoint wave functions in our construction. Furthermore, 
our approach is based on Painleve analysis while they obtain their DT by means of 
reductions. Noticing these differences, we are interested in constructing a DT for (|l]) by 
reducing the binary DT for a general linear problem. 
As in 0, we consider the following linear problem 

(d 2 + Fd + U)^ = tia 3 V, (24) 

where /i is the spectral parameter and 

F=( f » ^ 12 V U=( Uu Wl2 V a 3 =(l ° V (25) 
V hi /22 / \U 2 l u 22 ) \ -1 J 

If 

F = 0, u u = u 2 2 = u, U12 = u 2 i = v, (26) 

above linear problem is equivalent to the spatial part of the linear problem for the HS-KdV 
system given in last section. 
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The binary Darboux transformation has been constructed in general case, for example, 
contains a detailed discussion on this issue. In the present case, we consider 



L0 



8 y - a 3 (8 2 + Fd + U) 9 = 0, (27) 



where y is a new independent variable, F and U are those matrices given by It is 

easy to see that ( p4|) is a dimension reduction of (p7|) . 

What we will do next is to study reduction of the general binary DT and construct a 
proper DT for HS-KdV system (H). To construct a binary DT. we also need to consider 
the adjoint linear problem 

" - d y - (d 2 - OF 1 + E/*)<r 3 ] $ = 0. (28) 

Now we take as a matrix solution of the linear system ( U% and $ as a matrix solution 
of the adjoint linear system (|28|). Then we define a potential matrix by 

n x = $*e, 

here and in the sequel superscript * denotes matrix transposition. 
It is known that 

g = i - en -1 ^ 1 ^ 

satisfies the following equation 

{dy - a 3 (d 2 + F8+ U))G = G(d y - a 3 (8 2 + F8 + U)), 

where the hatted matrices F and U are the 2x2 matrices with the new transformed 
variables. Therefore, one has a binary DT (for wave functions) 

if = g^ 

in this general case. As argued in ||13|| , this general DT preserves certain properties of the 
operator. Therefore, we choose in the sequel the adjoint wave function as 



$ = —AO, A 



-1 

1 



This reduction is consistent with the DT. Furthermore, a simple calculation shows that 
under this reduction F is not alerted 

P - f = eir 1 ^ - crgeir^Vg = o, 

and the potential matrix Q is found to satisfy 

tt x = det(6)A (29) 



We now have a DT for the linear problem fl27Q with F — 0. Since the HS-KdV system 
(0) is a 1 + 1-dimensional system, we have to make further reduction. That is, we need to 
do dimensional reduction. Furthermore, U has to be symmetric and its diagonal entries 
have to be identical, namely, U need to take the form (|2T)|). To this end, we take 

I 02 01 J ' 
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where (61,62) is the solution of the following system with /j, — A 

\<h ) t \ v x -\u x J \ 2 J \ -2v u-2/i J \ 02 J x 

In this case to show that the general DT is reducible, one has to prove the following 
identity holds 

<j 3 (d 2 + u)(i - en^d^QA) = (1 - en-'d^QA)^ 2 + u), 

where U = ^ ^ ^ ^ . Above equation yields 

u-u = (eir 1 ^^ + ^(en^eA)^ + (bq-^qa - ^(eir 1 )^^, (30) 

and 



+ 



Gir^ArgG - asQQ^QAe^Q^d^QA = 0. (31) 



The validity of above identity fl3~T| ) can be checked by straightforward calculation. Thus, 
we have a DT for the linear problems (p4[): 

£ = (1 - QQ^d^OA)^, (32) 

and the transformation for U is given by (|30|). This shows that the general Darboux 
transformation can be reduced to the special case: the HS-KdV case. Of course, we also 
have to work out the time evolution part of the DT. This can be easily done by means of 
MAPLE. 

To compare this DT with the DT obtained in last section, we rewrite above results in 
terms of scalars rather than vectors. From the transformation for U (|30|), we readily have 

cj / 1 \ ~ r, 0l&2x - 6l x 6 2 

u - u = 2(lnp) xx , v-v = 2 , 

P 

where 

Px = 6\ - 6\, 

is derived from (|29[). The transformation (1321) for the wave functions is 



P 1 

h + ~, 02 = 02 + -, 
P P 



where 



62 

P = T~ (#201z - 6*2x01 + #1 4>2x ~ #1x02) + 

A + ix 

Q 

q = T ' (0201x - 02x01 + 01 02x - 01x02) + 

A + /i 



01 (a a. 


A- 


V^lS^lx 

p 


62 




A- 


-(010!, 
P 



For the temporal part, we need time evolution for p, which is 

pt = 2\{el + el) + \u(Q\ - e\) + e\ x - q\ x . 

With all these in hand, we can show that the transformed variables and wave functions 
solve the t— part of the linear problem. 

Therefore, we do find a DT for the HS-KdV system. The DT presented in this section 
seems different from the one obtained in last section, but indeed they are same. To 
establish the equivalence, we need to do the following transformations 

6l + 9 2 = fl, 01 ~ 02 = /2, U2 = U, V 1 = V, 

then an easy calculation shows that the DT obtained by reduction is nothing but the one 
constructed by Painleve analysis. 



4 Iterated DT and Exact Solutions 



In general, a DT can be iterated and a compact representation of solutions can be derived 
in terms of some special determinants such as Wronskian, Grammian. We will show that 
it is also the case for our DT. 

We notice that our DT for the wave functions flT3T) can be written as 



f]i4> 2 dx. 



Vi f x ~ V2 

01 = 01 / 4>lV2dx, 02 = 02 

T J T 

To do iteration, we take iV solutions of the linear system (f?|-|T0|) {rji \ 77J5 ) with respect to 
Xk (k = 1, . . . , N). After iterations, we have the new wave functions 

J x r)i^ 772*' dx ■ ■ ■ J x i]^ 7/2^ dx ry} 1 ' 



1 



-[N] 



J x </)irj [ 2 ] dx 



J x rji^rj^dx r/i^ 
I x 4>iV [ 2 N] dx 0i 



(33) 



and 
where 



r[N] 



rv?v [ 2 1] dx 



I x Vi N] V2 ] dx 



V2 



F rj?T] [ 2 N] dx rvPfodx 



J x r]^i] 2 N ^dx J x i]f l (p2dx 

(p 2 



V [ 2 N] 



u[N] =u + 2(\nr[N]), J 



t[N] = det 



[fc] [i] 

m m 



(34) 



Although after a single DT, the field variable v has a very nice representation, we do not 
have a similar one after iterations. However, if we put u = w x , we may solve the first 
equation of (P and obtain 



6 



[w xxx + 310^ - 2w t 
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In this way, we can find the transformation for the field v. Therefore, the solutions for 
(0) are represented in terms of a Grammian r[N}. 

As a final part of this section, we generate solutions for the HS-KdV system ([!]) by 
means of our DT. Let consider the simplest case: u 2 = 0, V\ = c (c is an arbitrary 
constant). Thus we need to solve 

Vixx = -crji - i\rj 2 , rj 2xx = cr] 2 + iXrjt, 
rj lt = -2cr) lx - 2i\r) 2x , r) 2t = 2cr) 2x + 2i\r) lx , 

where for convenience we made a scaling for the spectral parameter: A — > iX (i 2 = — 1). 
The general solutions of the above system read as 



Ci 



e ^/n(2Kt+x) _|_ c ^ e -^K(2nt+x) _|_ c ^ e -i^K(-2nt+x) _|_ c ^ e i^K(-2Kt+x) 



V2 = Cl (c + K)e^ 2Kt+x) + c 2 (c + n)e-^ 2Kt+x) 

+ C 3 (C - K ) e -^{-^t+x) + _ K y^{-2Kt+x) 

where k 2 = c 2 + A 2 , Cj (j — 1, ... ,4) are arbitrary constants. Solving 



(35) 
(36) 



r x = V1V2, n = 2i\(rjl - r?f) - 2r] lx r) 2:i 



we find the potential 

i 

T 



\/ K 2 — C 2 



24H 



12k [c\c 2 (c + k) + c 3 C/i(k — c)] t + 2 [cic 2 (c + k) — c 3 c 4 (/t — c)] X 

+ 



C 2 (c + k) c 2y^(2 K f+a) _ C 2 (c + k) c - 2y /j:( 2Kt+x ) Zc|(c ft) ^2i^,{-2Kt+x) 



2^R 



2^ 



2v^ 



|_ e 2j % /^(-2Kt+a;) _ (1 + ^JC2C4C e - v ^[2(l+i)/tt+(l-i) a ;] 



,1 z)c 2 C3C c _^ [2(1 _. )Kt+(1+ . ):r1 (1 + ^)CiC 3 C cx/ ^ [2(1+ . )Kf+(1 _. )x1 



K 



K 



i)x] 



(37) 



and the solutions for (|l|) are 

u = (lnr)^, 



v — c + 



f]\xf]2 - VlV2x 



where rji, r\ 2 and r are given by (|35"D, (|36|) and (|37D respectively. Our solutions contain 
the known 1-soliton solutions as special cases. Indeed, by setting c\ = c 3 = 0, we arrive 
at 



ic 2 2 



(C+K) 



-2^/~K,{2Kt+x) 



V1XV2 ~ VlV2i 



2 \ k(k — c) 

fjj^ -2^(2rt +I )i 



1 2(1 + ?)cC 4 



c 2 (c + Kj 



2C2 C — M 



.2(9 



Cf(c+ ft) 



\/ K 2 — C 2 



where # = -y/fi;[2(l — + (1 + if we further put c = 0, we obtain the 1-soliton 
solution first found by Hirota and Satsuma M. We remark that our DT allows us to 
construct multi-soliton solutions for HS-KdV system ([!]). 
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5 Conclusions 

In this paper, we constructed a new DT for the HS-KdV system. We notice that our 
DT is a binary DT while the one constructed by Leble and Ustinov is a classical DT. 
The situation is similar to the celebrated KdV case. We obtained our DT by using two 
methods, namely Painleve analysis and reduction method. In general, one may obtain a 
Backlund transformation if the wave functions involved in a DT are eliminated. Thus, 
one may ask if it is possible to obtain a Backlund transformation from our new DT. We 
realize that it is not trivial to do this sort of elimination in the present case. On the other 
hand, we remark that DT is easier to handle than Backlund transformation in order to 
construct explicit solutions. 
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